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I. INTRODUCTION 

The question whether a given state cannot be de¬ 
scribed within a classical theory has been considered 
one of the most serious since the early days of quantum 
physics [TH3] (for a review see, e.g., Ref. 0D- Nonclassi¬ 
cality and entanglement, which is one of the nonclassical¬ 
ity manifestations, are the most important properties of 
optical fields studied in quantum optics. Such fields have 
no classical analogs and as such they have been found 
interesting for many reasons. Nonclassical properties of 
such fields have been found useful both for elucidating the 
principles of quantum mechanics and in various applica¬ 
tions including, e.g., quantum information processing [5], 
quantum metrology Era, and highly-sensitive measure¬ 
ments 0. 

From both the theoretical and the experimental points 
of view, the nonlinear process of parametric down- 
conversion, in which photon pairs are generated, has 
played an important role here from the beginning of in¬ 
vestigations gUME]. Its individual photon pairs have 
been exploited in many fundamental experiments testing 
nonclassical behavior predicted by quantum physics H 
HU. It has also allowed the generation of more intense 
fields having their electric-field amplitude quadratures 
squeezed below the vacuum level HSHH3, exhibiting sub- 
shot-noise correlations Hi OS or having sub-Poissonian 
photon-number statistics [ 201422 ] . 

In quantum optics, the definition of nonclassicality 
is based upon the Glauber-Sudarshan P representa¬ 
tion mmm of the statistical operator of a given 
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field. The commonly accepted formal criterion for dis¬ 
tinguishing nonclassical states from classical ones is ex¬ 
pressed as follows pi Hunan)]: A quantum state is 
nonclassical if its Glauber-Sudarshan P function fails to 
have the properties of a probability density. Alterna¬ 
tively, several operational criteria for nonclassicality of 
either single-mode pHMTl or multimode [2H14T?(7j fields 
have been revealed. Their derivations are based either on 
field’s moments [Ml 130 SB or on direct reconstruction of 
quasidistributions of integrated intensities [52HM] . Also, 
criteria derived from the majorization theory have been 
found [35( [36]. 

Entanglement (or inseparability) is a special nonclassi¬ 
cal property that describes quantum correlations among 
(in general) several subsystems that cannot be treated by 
the means of classical statistical theory m- Various ap¬ 
proaches have been developed for discrete and continuous 
variables to reveal entanglement. This property has been 
exploited in suggesting an entanglement criterion and the 
related entanglement measure (referred to as the negativ¬ 
ity) based upon the partial transposition of a statistical 
operator 1351441 ] . Another approach has been based on 
the violation of the Bell inequalities written for different 
mean values including the measurement on both parts of 
a bipartite system [Ml- Also, a method using positive 
semi-definite matrices of fields’ moments of different or¬ 
ders 23 033 bas been found to be very powerful. We 
would like to stress at this point that entanglement is 
a very crucial tool in today’s quantum information pro¬ 
cessing. 

In this contribution, we study nonclassicality by apply¬ 
ing the Lee nonclassical depth [45] as well as entangle¬ 
ment via the negativity maim for (in general) noisy twin 
beams of different intensities. Such fields occur under real 
experimental conditions in which a nonlinear crystal gen¬ 
erates both photon pairs and individual single photons 
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(noise). Nevertheless, the signal and idler fields together 
form a bipartite quantum system. We note that entan¬ 
glement and nonclassicality of twin beams generated by 
down-conversion seeded by thermal light have been ana¬ 
lyzed in Refs. (TTH - TTSl . In this case, noise present in the 
incident thermal fields participates in the nonlinear pro¬ 
cess and generation of photon pairs. This weakens its 
detrimental effect on entanglement and nonclassicality of 
twin beams and allows us to have entangled twin beams 
with a larger amount of noise. 

Here we also study the problem of entanglement di¬ 
mension via the negativity N for general twin beams and 
the Schmidt number I\ for noiseless twin beams in a pure 
state. Namely, we estimate how many degrees of free¬ 
dom of two fields comprising a twin beam are entangled 
based on the results of Ref. [49] for axisymmetric states. 
On the other hand, the participation ratio R s [5U1 de¬ 
termined from the reduced statistical operator p s of the 
signal (or idler) field gives the number of degrees of free¬ 
dom in this field serving to describe both entanglement 
and noise. It varies from R s = 1 (for a pure state p s ) 
to R s = d = dim(/5 s ) for the completely mixed state 
p s = I/d. We note that the participation ratio R s gives 
an effective number of states in the mixture p s implied 
by the property that it is a lower bound for the rank of 
p s . Moreover, the logarithm of R is the von Neumann- 
Renyi entropy of second order |50| . The inverse of the 
participation ratio is referred to as the purity (or lin¬ 
ear entropy). Various methods for direct measuring the 
Schmidt number K (even without recourse to quantum 
tomography) were proposed for noiseless twin beams (see, 
e.g., Refs. [SMS)- The method of Ref. [53] was recently 
realized experimentally [56] . We note that the negativity 
can also be measured without applying quantum tomog¬ 
raphy as described, e.g., for two polarization qubits using 
linear optical setups mm- 

The paper is organized as follows. In Sec. II, the model 
of parametric down-conversion providing an appropriate 
statistical operator of a twin beam is presented. Entan¬ 
glement of the twin beam is addressed in Sec. Ill using 
the negativity. The nonclassical depth is introduced in 
Sec. IV to quantify nonclassicality. The relation between 
the negativity and the nonclassical depth is also discussed 
in Sec. IV. The dimensionality of a twin beam described 
by the participation ratio together with the entanglement 
dimensionality described by the negativity is analyzed in 
Sec. V. Properties of M-mode twin beams are discussed 
in Sec. VI. Section VII is devoted to experimental multi- 
mode twin beams containing also noise embedded in in¬ 
dependent spatiotemporal modes. Conclusions are drawn 
in Sec. VIII. 


II. QUANTUM MODEL OF A TWIN BEAM 

To describe the generation of a single-mode twin beam 
by parametric down-conversion, we adopt the approach 
based on the Heisenberg equations derived from the ap¬ 


propriate nonlinear Hamiltonian H llA m, 

Hint = -hgAiA 2 exp(iwt - i0) + H.c., (1) 

where Ai (A{) and A 2 (Hj) represent the annihilation 
(creation) operators of the signal and idler field, respec¬ 
tively, and g is a real coupling constant that is linearly 
proportional both to the quadratic susceptibility of a 
nonlinear medium and to the real pump-field amplitude. 
The interaction time is denoted t , w ( cf >) is the pump-field 
frequency (phase), and ui\ and ui 2 stand for the signal- 
and idler-field frequencies, respectively. The law of en¬ 
ergy conservation provides the relation uj = uii+lu 2 . H.c. 
is the Hermitian conjugated term. In a real nonlinear 
process, also noise occurs. It can be described by the 
Langevin forces L belonging to a reservoir of chaotic os¬ 
cillators with mean number of noise photons (n d ). 

The Heisenberg-Langevin equations corresponding to 
the Hamiltonian H ult are written as 

— (iuq + 71 )Ai + i gA\ exp(—ic ut + i</>) + Li, 

— (i oj 2 + 72 ) 4)2 + i gA\ exp(—i ut + i 4>) + L 2 , 

( 2 ) 

where the constant 71 (72) describes damping in the sig¬ 
nal (idler) field. The Langevin operators Li (for i = 1,2) 
have the properties 

(L i ) = (L\)=0,(LlL j )=2 lj (n d )S lJ , 

(Lit]) = 27 j ((n d ) + 1) 5ij, (3) 

where Sjj stands for the Kronecker symbol. 

Using the interaction representation [Aj (t) = 

cij(t) exp(— icojt)] and neglecting damping together with 
the Langevin forces, the solution of Eq. ([ 2 ]) attains the 
form 


d A! 
d t 

di 2 


ai(f) = di(0)u(t) + id\(0)v(t) exp(i(j)), 
a 2 (t) = a 2 (0)u(t) + ia{(0)v(t) exp(i0), (4) 

in which u(t) = cosh (gt) and v(t ) = sinh(gt). 

Statistical properties of the twin beam are then de¬ 
scribed by the normal characteristic function C[v defined 
as 


CV(/31,/3 2 ) = Tr 


/3exp(/3ia{ + ^ 2 ^ 2 ) exp(—/3j(ai 


/ %a 2 ) , 


(5) 

where Tr denotes the trace. Using the solution given in 
Eq. Q, the normal characteristic function CV attains 
the Gaussian form [591 . 


CV(/3i>^ 2 ) — ex P [— {\(3i\ 2 Bi + \/3 2 \ 2 B 2 ) + Di 2 f3l^2 + 
+■^12^1^2] > ( 6 ) 


in which /3i and (3 2 denote independent complex vari¬ 
ables. For the undamped and noiseless case, we have 
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D 12 = (Aai A 02 ). Also the mean number B p of the gen¬ 
erated photon pairs is determined as B p = (Aa{ A di) = 
(AaJ A a 2 ). When damping and noise are also consid¬ 
ered [52], the parameters B a (for a = 1,2) contain addi¬ 
tional noise contributions characterized by the parame¬ 
ters B s and B{, i.e., Bi = B p + B s and B 2 = B p + B[. 
Whereas the parameter B p gives the mean number of 
photon pairs, the parameters B s and Bi correspond to 
the mean number of noise photons coming from the 
signal- and idler-field reservoirs, respectively. On the 
other hand, the parameter D 12 describing mutual cor¬ 
relations between the signal and idler fields is not influ¬ 
enced by the noise since |Z?i 2 1 2 = B p (B p + 1). 

The statistical operator p of the twin beam then ac¬ 
quires the form El 


P = \ 2 J d 2 fad 2 (3 2 Ca(Pi, P 2 ) ■ exp (^2 ~ j : 

(7) 

In Eq. §, C A (P u p 2 ) = Cj/(Pi,p 2 ) exp(—|/?i | 2 - |/? 2 | 2 ) 
denotes an anti-normal characteristic function and sym¬ 
bol : : means normal ordering of field operators. 

Performing integration in Eq. ([T]) we express the sta¬ 
tistical operator p in the form 


P = 


K 


: exp 


B 2 „ B !„■(■., 

— —a\ a± -— a.,a 2 

K K 



\B> 12 /G . , 

- 1 aia 2 + a[a' 2 I 

I\ ' ' 

■ 1 

( 8 ) 

12 2 . The parameters B a 

intro- 


where K = B±B 2 — \Di 
duced in Eq. Q are related to anti-normal ordering of 
field operators and are given as B a = B a +1 with a = 1,2. 
Decomposing the statistical operator p in the Fock-state 
basis we finally arrive at the formula 


Pij,kl 


00 n p r 

EEEEH) 


n—0 p— 0r—0 t—0 


n—r 


B2- p B'r r k- n - 1 

(n — p)\(p — r)l 


x 


\D12Y 

(r — 1)\t\ 


( ij\a 


]n-p+t ~]p-r+t ~n-p+r-t ~p- 

1 a 2 a l a 2 


‘l kl). 


(9) 


Direct inspection of Eq. for the matrix elements of 
the statistical operator p written in Eq. reveals that 
all nonzero elements can be parameterized by only three 
indices, 


Pi,j,i-\-d,j-\-d — 


1 /(i + d)\(j + d)\ 


ma x(i,j) 


K 

to ! 


j! 


E c ,i 


m =0 


(to + d)\ 


f\D 12 \ 


r 1 -! 


V k 


d+2m 


, (io) 


assuming d > 0. Moreover, pi^ i+d j +d = pi+d,j+d,i,j, 
X a = 1 — B a /K with a = 1, 2, and C l m and C denote 
the binomial coefficients. 


III. NEGATIVITY OF THE TWIN BEAM 


The negativity TV of a mixed bipartite system defined 
on the basis of the Peres-Horodecki criterion for a par¬ 
tially transposed statistical operator [38] ISS E] is useful 
for quantifying the entanglement of the twin beam. It 
can be expressed as 


N(p) = 


lld r ||i-l 


( 11 ) 


using the trace norm ||p r ||i of the partially transposed 
statistical operator p T . The negativity essentially mea¬ 
sures the degree at which p T fails to be positive. As such 
it can be regarded as a quantitative version of the Peres- 
Horodecki criterion for separability [38][39]. According to 
Eq. (11), the negativity TV is given as the absolute value 
of the sum of the negative eigenvalues of p T . It vanishes 
for separable states. It is worth noting that the negativity 
TV is an entanglement monotone and so it can be used to 
quantify the degree of entanglement in bipartite systems. 
Moreover, the negativity does not reveal bound entangle¬ 
ment (i.e., nondistillable entanglement) in systems more 
complicated than two qubits or qubit-qutrit m- 

To determine the negativity TV we consider the eigen¬ 
value problem for the partially transposed statistical op¬ 
erator p T . The statistical operator p T expressed in the 
Fock-state basis attains a characteristic block structure. 
The smallest block has dimension 2 and each successive 
block has dimension larger by 1. For a given M one has 
a block of dimension M + 1. Such a block represents a 
matrix of M + 1 isolated states; the sum of indices of 
their statistical operators equals 2 M, 


P0M,0M Po M—1,1 M ■■■ Poo,MM 


~r _ PlM,0M-l Pi M—1,1 M—1 

Pm — 

Pm m, 00 


Pmo,mo 


( 12 ) 

It can be shown that eigenvalues of a block 
of dimension M + 1 can be expressed as 


i/_,..., V 5 using the eigenvalues 


'+ > + 


zz+ and V- of a block with dimension 2 : 


v± = 1-=- 

2 K 


Bi + B 2 T \ [B 2 - B 


2 


4| D 


121 


, (13) 

The negative eigenvalues can only be those containing 
odd powers of They form a geometric progression 
whose elements can be summed to arrive at the formula 
for the negativity TV: 


TV = 


1 3(Hi + B 2 ) + J(B 1 - B 2 y + 4|D 12 | 2 - AK - 2 


AK - 2(1?! +B 2 ) + 1 


CM) 


Expressing parameters Bi, B 2 , and |Di 2 | 2 in Eq. (14) 


in terms of parameters B p , B s , and B u we arrive at the 
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FIG. 1: (Color online) Negativity IV as a function of the mean 
photon-pair number B p for noiseless twin beams (i.e., B s = 
B\ = 0)) according to Eq. (161. 


formula 

N = {2B p -(B s + B i )(^B p + l)-4B s B i 

+ yj(B b - A) 2 + AB P (B P + 1)} 

x {4 (B a + Si)( 2 B p + 1) + 8 B s Bi + 2}” 1 .(15) 

Equation © simplifies considerably for noiseless twin 
beams: 


N = B p + B p (B p + 1). 


(16) 


According to Eq. (16), all noiseless twin beams are 


entangled. The more intense the noiseless twin beams 
are, the more entangled the signal and idler fields are 
(see Fig. [I]). The presence of noise in a twin beam can 
even completely destroy entanglement, as the analysis of 
Eq. © shows. Indeed, the condition N > 0 for entan¬ 
glement can be rewritten using Eq. (15) as follows: 


B p [l-(B B + Bi)] >B S B { . (17) 

Condition © cannot be fulfilled for any value of B p 
provided that B s + B\ > 1 . Thus, the twin beam can be 
entangled only when 


B s + B { < 1 


and 


B p > 


B s Bj 

1 ~(B S + Bi )' 


(18) 


The behavior of the negativity N of noisy twin beams 
dependent on the noise parameters B s and Bi is illus¬ 
trated in Fig. [2] for several values of the mean photon- 
pair number B p . It holds in general that the greater the 
value of the mean photon-pair number B p , the greater 
the value of the negativity N. This can be explained as 
follows. The more intense twin beams, with their ther¬ 
mal statistics, are effectively spread over a larger number 



FIG. 2: (Color online) Negativity N, given in Eq. (15), as 
a function of the mean noise photon numbers B s and Bi in 
the signal and idler modes, respectively, assuming the mean 
photon-pair number B p equal to 0.5 [bottom light-gray (yel¬ 
low) area], 1 [gray (green) area], 2 [dark-gray (blue) area] and 
4 [top, black area]. The larger B p , the larger the negativity 
N. 


of the Fock states. This naturally results in the larger ef¬ 
fectively populated Hilbert spaces used to describe the 
entanglement. The greater value of the negativity N 
means a greater effective number of the paired modes 
building the entanglement, i.e., a greater value of the en¬ 
tanglement dimensionality, as defined in Sec. V. Also, the 
greater the value of the mean photon-pair number B pi 
the larger the amount of overall noise B s + B[ acceptable 
in an entangled twin beam (see Fig. [3]). The curves plot¬ 
ted in Fig. [^indicate that entanglement is more resistent 
to noise when the noise is distributed in the signal and 
idler fields asymmetrically. We note that separable states 
(i.e., with IV = 0) contain, in general, paired, signal, and 
idler noisy contributions. However, the noisy contribu¬ 
tions are sufficiently strong to suppress the “entangling 
power” of the photon-pair contribution and so the state 
effectively behaves as a classical statistical mixture of the 
signal and idler fields. 

The decomposition of the partially transposed statis¬ 
tical operator p r into blocks in its matrix representation 
and the fact that a block (subspace) with dimension M+ 1 
describes only states with up to M photons in the signal 
(and also idler) field can be used to define the distribu¬ 
tion d/v of the negativity N fulfilling the normalization 
condition 

OO 

d N (M) = N. (19) 

M—l 

For a given M, the element djsr(M) of this distribution 
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FIG. 3: (Color online) Curves giving the boundaries between 
entangled and separable twin beams and determined accord¬ 
ing to Eq. (181 plotted in the plane spanned by the mean noise 
photon numbers B s and Bi assuming the mean photon-pair 
number B p equal to 0.01 [dotted (red) curve], 0.1 [dash-dotted 
(yellow) curve], 0.5 [dashed (green) curve], 2 [long-dashed 
(blue) curve], and B p = 100 [solid black curve]. Entangled 
states are localized in the lower-left corner of the plane. The 
larger B p , the larger the area containing entangled states. 


is given as the sum of the absolute values of the negative 
eigenvalues belonging to the block of dimension M + 1. 
The distribution dw of the negativity provides insight 
into the internal structure of entanglement. It tells us 
how entanglement is distributed in the Liouville space 
of statistical operators. Typical distributions g^v of the 
negativity for noiseless as well as noisy twin beams are 
plotted in Fig. [4] A teeth-like structure occurs for smaller 
numbers M in noiseless twin beams. Noise tends to sup¬ 
press this structure, as is evident from the comparison 
of the distributions d n plotted in Figs. |4][a) and[4][b). 
We note that the densities of the negativity have already 
been introduced for bipartite entangled states composed 
of a qubit and continuum of states [60] [BT] as well as two 
continua of states. 


IV. NONCLASSICAL DEPTH OF THE TWIN 
BEAM 

To quantify nonclassicality of the twin beam we apply 
the nonclassical depth r [JS] derived from the thresh¬ 
old value Sth of the ordering parameter at which the 
joint signal-idler quasidistribution of integrated intensi¬ 
ties becomes nonnegative [2H UjI] . We adopt the defini¬ 
tion r = (1 — s t h)/2. We note that the joint signal-idler 
quasidistribution of integrated intensities attains nega¬ 
tive values for 1 > s > Sth for which r > 0. The 
threshold value s t h can easily be obtained from the con¬ 
dition ([A(VF S — W\)} 2 ) = 0, which determines the point 



FIG. 4: (Color online) Distribution djv of negativity N given 
in Eq. (191 assuming B p = 2 and (a) B a = Bi = 0 and (b) 
B s = Bi = 0.1. Note that —djv(M) corresponds to the sum of 
all the negative eigenvalues for the (M + 1)-dimensional block 
of the partially transposed statistical operator p r . Thus, 
disr(M) shows the internal structure of entanglement in the 
Liouville space. 


of the transition between quantum and classical single¬ 
mode twin beams [34] , This results in the following for¬ 
mula for the nonclassical depth r: 


1 

T = - 
2 


B s - Bi) 2 + AB p (B p + 1) - 2 B p - B s - B { ) 


( 20 ) 


Assuming noiseless twin beams, Eq. (20) simplifies to 


r — \ B p (B p + 1) — B p . 


( 21 ) 


According to Eq. (21), all noiseless twin beams are non¬ 


classical. The greater the mean photon-pair number B p , 
the greater the value of the nonclassical depth r (see 
Fig# This depth r reaches its greatest value, 1/2, in 
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the limit of an infinitely intense twin beam ( B p —► oo). 
We note that r = 1/2 corresponds to symmetrical order¬ 
ing of the field operators. 


On the other hand, and according to Eq. (201, noise 


only degrades nonclassical behavior of a twin beam, as 
documented in Fig. [gJ If the noise is equally distributed 
in the signal and idler fields ( B s = B^), the nonclassical 
depth r determined along Eq. ( [20| gives the mean num¬ 
ber B s + B[ of noise photons needed for suppressing the 
nonclassicality of the twin beam. So, the larger the value 
of the nonclassical depth r is, the more nonclassical the 
field is. On the other hand, formal application of Eq. (20) 


to classical noisy twin beams results in negative values of 
the nonclassical depth r. Their absolute value |r| can be 
considered a measure of classicality of noisy twin beams 
in the sense that it quantifies the mean number of photon 
pairs needed to transform a classical twin beam into the 
classical-quantum boundary r = 0. 

Condition t = 0 for the transition from quantum to 
classical twin beams applied to Eq. (20) results in the 


same relation among parameters B p , B s , and B\ as de¬ 
rived in Eq. © for the boundary between entangled and 
separable twin beams. Thus, entangled twin beams are 
nonclassical, whereas separable twin beams are classical. 
This means that nonclassical twin beams may contain on 
average only less than one noise photon ( B s + B[ < 1). 
We note that inequality ( p~7|) represents the Simon cri¬ 
terion for nonclassicality of Gaussian states as shown in 
Ref. (62]. 


Comparison of Eqs. (16) and (21) made for noiseless 


twin beams reveals a simple relation between the nega¬ 
tivity N and the nonclassical depth r: 


N = 


1-2 T 


( 22 ) 


Direct calculation based on Eqs. (15) and (20) then con 


FIG. 6: (Color online) Nonclassical depth r given in Eq. ( |20| ) 
as a function of the mean noise photon numbers B s and Bi 
for the mean photon-pair number B p equal to 0.1 [bottom, 
light-gray (yellow) area], 0.5 [gray (green) area], 4 [top, dark- 
gray (blue) area]. The greater the value of B p , the greater 
the value of r. 


firms that relation (221 holds even for a general noisy 
twin beam. We thus have a one-to-one correspondence 
between the value of the negativity N and the value of 
the nonclassical depth r. Moreover, according to Eq. (22) 
the negativity N is an increasing function of the nonclas¬ 
sical depth r, and vice versa (see Fig. [7]). There exists a 
deep physical reason for this correspondence. The non¬ 
linear process emits photons in pairs into the signal and 
idler fields, which creates entanglement between these 
fields. It is this entanglement that gives rise to nonclas¬ 
sical properties of twin beams, as the classical statistical 
optics is unable to describe pairing of photons appropri¬ 
ately. 


V. DIMENSIONALITY OF THE TWIN BEAM 

Three different numbers are needed to determine the 
dimensionality of a general noisy twin beam. The di¬ 
mensionality Kent of entanglement gives the number of 
degrees of freedom constituting the entangled (paired) 
part of the twin beam. We also need additional degrees of 
freedom to characterize the noisy parts of the twin beam. 
As the amount of noise is, in general, different in the sig¬ 
nal and idler fields, we have independent participation 
ratios R s and i?i for both fields. The entanglement di¬ 
mensionality Kent for bipartite states with axisymmetric 
statistical operators can be given in terms of the nega- 
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FIG. 7: (Color online) Negativity N as a function of the non- 
classical depth t, according to Eq. (]22[). 


tivity TV by a simple formula : 

K ent (p) = 2N(p) + l = \\p r \\ 1 . 


(23) 


Strictly speaking, it is the least integer > I \ ent that gives 
a lower bound to the number of entangled dimensions 
between entangled subsystems (paired modes) of p [49] . 
According to Eq. (23), the entanglement dimensionality 


K ent equals 1 for separable states (TV = 0). It linearly 
increases with the negativity N. As the noise described 
by the mean noise photon numbers B s and B, decreases 
the values of the negativity N, it also decreases the val¬ 
ues of the entanglement dimensionality K en t- We note 
that, for pure states, the Schmidt number is also a good 
quantifier of the entanglement dimension K ent |63H65| . 
The Schmidt decomposition of pure states accompanied 
by convex optimization can even be applied for quan¬ 
tifying the entanglement dimension of mixed entangled 

states [37]. 

On the other hand, the noise present in the signal and 
idler fields requires additional degrees of freedom for its 
description. These degrees of freedom are, together with 
those reserved for describing entanglement, determined 
by the participation ratios R s and R\ derived from the 
signal- and idler-field reduced statistical operators p s and 
Pi, respectively [64, 66): 


Rn = 


1 


44’a \Pa 


a = s, i. 


(24) 


Using Eq. (10) the matrix elements p s jj can be written 
as 


Ps,jj — 



(26) 


Substituting Eq. (26) into Eq. (25) we obtain a simple 


formula for the participation ratio R s \ 
Rs = 2(-B p + B s ) + 1. 


(27) 


The same considerations made for the signal field apply 
also to the idler field. 

To find the relation between the entanglement dimen¬ 
sionality .Kent and the participation ratios R s and R\ we 
consider for a while the noiseless twin beams in pure 
states. In this case, the elements p s _ r j of the reduced 
statistical operator p s , written in Eq. (26), immediately 
give the squared Schmidt coefficients m Combining 
Eqs. (Jl6|, ([23]) , and ([27]) we arrive at the formula 


K ent = Rs + VRs - 1- 


(28) 


Equation (28) shows that, excluding weak noiseless twin 
beams, A" ent ~ 2 R s . This means that the definitions of 
the entanglement dimensionality and participation ratio 
set different boundaries for the Schmidt coefficients Cj 
included in the approximative description of a noiseless 
twin beam with the wave function 


= c ib%li)i- 

3=0 


(29) 


Using Eq. (10), the coefficients Cj in Eq. (29) are obtained 
in the form 


c,- = 


V 


(B p + 1) J+1 ’ 


(30) 


which is in agreement with the thermal photon-number 
statistics of the signal (or idler) field. We note that 
the ratio cx mt - i/cr b _i of boundary coefficients is given 
by the expression [B p /( 1 + B p )] B ? +1 . When B p —> oo 
CK eDt -i/c Rs -i -> 1/e. 

To compare the values of entanglement dimensional¬ 
ity and the participation ratio for general twin beams we 
have to eliminate the effect of different boundaries set 
by different definitions, as revealed by considering the 
pure states. Using the formulas derived for noiseless twin 
beams, we introduce the modified entanglement dimen¬ 
sionality Kent as follows: 


Equation (101, giving the matrix elements of the statis¬ 


tical operator p, guarantees a diagonal form of the re¬ 
duced statistical operators p s and p i of the signal and 
idler fields, respectively. In this case, Eq. (24) can be 
rewritten in the form 


R s = 


l 


E : p: 


,2 ' 
s,33 


(25) 


Kent. 


2 B n 


2B p + 1 


1 + 2 JBI 


K fil 


B n 


(31) 


Definition (31) of the modified entanglement dimension¬ 


ality Kent guarantees that the values of modified entan¬ 
glement dimensionality K ent and participation ratios R s 
and Ri of noiseless twin beams are equal. 





























FIG. 8: (Color online ) M odified entanglement dimensional¬ 
ity A' ent given in Eq. |3l| [lower, dark-gray (blue) area] and 
signal-field participation ratio R s given in Eq. (271 [upper, 
gray (red) area] as they depend on the mean noise photon 
numbers B s and B\ assuming the mean photon-pair number 

B P = 1. 



FIG. 9: (Color online) Coefficient r en t given in Eq. ( |32| ) ver¬ 
sus the mean noise photon numbers B a and B\ for the mean 
photon-pair number B p equal to 1 [upper, dark-gray (blue) 
area], and 10 [lower, gray (red) area]. 


The values of the modified dimensionality K ent of en¬ 
tanglement and the signal-held participation ratio R s are 
compared in Fig. [8] for the mean photon-pair number 
B p = 1. Whereas the values of the modified entangle¬ 
ment dimensionality K ent decrease with increasing val¬ 
ues of the mean noise photon numbers B s and A;, the 
values of the signal-held participation ratio R s increase 
with increasing values of the mean signal-held noise pho¬ 
ton number B s . We note that the values of the signal-held 
participation ratio R s are greater than those of the mod¬ 
ified entanglement dimensionality K ent even for B s = 0, 
as the presence of noise in the idler held (A; > 0) de¬ 
grades entanglement. 

The relative contribution of the degrees of freedom 
used for describing entanglement in a twin beam is an 
important characteristic. This contribution can be quan¬ 
tified via the coefficient r ent dehned as follows: 


T ent 


2A ent 
R s + i?i 


(32) 


the larger dimensions of the effectively populated Hilbert 
spaces of more intense twin beams and, thus, the more 
complex structures of their entanglement. As a conse¬ 
quence, relatively higher numbers of degrees of freedom 
serving to describe entanglement in more intense noise¬ 
less twin beams are “released” by the noise and enlarge 
the noise parts of twin beams. 

Alternatively to the participation ratio R, we may ap¬ 
ply the von Neumann entropy S' of a reduced statistical 
operator. Taking into account the diagonal form of the 
signal-held reduced statistical operator p s with the ele¬ 
ments written in Eq. (26), the signal-held entropy S s is 


in general determined along the formula 

S s = -Tr(/5 S In p s ) = - ^ p Stjj In (p s ,jj). 


(33) 


Considering the specihc form of matrix elements p s jj 
given in Eq. (26), the formula for entropy S s attains the 
form 


As shown in Fig. |9j the greater the values of the mean 
noise photon numbers B s and B \, the smaller the val¬ 
ues of the coefficient r en t- The comparison of surfaces 
of the coefficient r en t drawn for the mean photon-pair 
numbers B p = 1 and B p = 10 in Fig. [9] reveals seem¬ 
ingly paradoxical behavior. The values of the coefficient 
r en t decrease with increasing values of the mean photon- 
pair number B p . This behavior, however, naturally orig¬ 
inates in fragility of entanglement with respect to the 
noise. More intense twin beams (with greater values of 
B p ) are less resistant to a given amount of noise com¬ 
pared to low-intensity twin beams. This is explained by 


S s — (1 + B p + B s ) ln(l + Bp ■+■ B s ) 

-(A p + A s )ln(A p + A s ); (34) 


In stands for natural logarithm. Combining Eqs. (25) and 
(34), the entropy S s is revealed as an increasing function 
of the participation ratio R s : 


S s =l [(Rs + 1) ln(i? s + 1) - (R s - 1) In (As - 1)] - 1. 

(35) 

Analogous formulas for the idler-field entropy S) can eas¬ 
ily be derived. The general dependence of entropy S s on 
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the participation ratio R s is plotted in Fig. [To} We would 
like to note that the entropy S serves as a good measure 
of the entanglement for pure states. 


VI. TWIN BEAM COMPOSED OF M MODES 


In real experiments, twin beams are rarely composed 
of only one paired spatiotemporal mode H3D] [33]. We 
note that a twin beam composed of one paired mode 
represents an ideal field from the experimental point of 
view EH. For this reason, we consider a multimode 
twin beam containing M independent identical single¬ 
mode twin beams. Its statistical operator pM is given 
as Pm = ®mP using the statistical operator p written 
in Eq. Q. There are four parameters characterizing 
the twin beam: number M of modes, mean photon-pair 
number B p , mean signal-field noise photon number B s , 
and mean idler-field noise photon number B\. We note 
that such an M-mode twin beam represents a good ap¬ 
proximation of a real twin beam when all spatiotemporal 
modes participating in the nonlinear interaction are de¬ 
tected. 

The considered physical quantities behave differently 
with respect to the number M of modes. It has been 
shown in Refs. [59] and [34] that the nonclassical depth 
r does not depend on the number M of modes. On 
the other hand, the multimode negativity Nm, Nm = 
(1 + 2 N) m , as well as the participation ratios Rm,o., 
Rm,<i = Ra 1 f° r a = s, i, are multiplicative. We note 
that the form of the multimode negativity originates in 
the multiplicative property of the trace norm and its rela¬ 
tion to the negativity expressed in Eq. © m- In fact, 
the multimode negativity Nm coincides with the entan¬ 
glement dimensionality K ent defined in Eq. (231 for a 


single-mode twin beam. The multimode entropies Sm,o, 
a = s, i, are then additive. To reveal similar relations 
among the studied quantities as has been done for single- 



FIG. 10: (Color online) von Neumann entropy S's as a function 
of the participation ratio R a according to Eq. <|35|). 


mode twin beams, we have to define suitable quantities 
derived from those considered above. Defining the log¬ 
arithmic negativity N^ s = In (Nm) and the logarithmic 
participation ratios Rm 0 = ln(I?M,a)i a = s, i, we replace 
the multiplicative quantities with the additive ones. In¬ 
troducing the logarithmic negativity A/", logarithmic par¬ 
ticipation ratios 1Z 1 £ S , and entropies S a related per one 
mode, 


AT = 

Ra = 
S„. = 


N 


log 

M 


M 


R 


,log 

M,a 


M 

R]\4,a 

M 


= ln(l + 2 V), 

= HRa), 

= s a , 


(36) 


with a = s, i, we reveal the suitable quantities. The quan¬ 
tities defined in Eq. (36) together with the nonclassical 
depth r behave qualitatively in the same way as those 
defined for single-mode twin beams discussed above. Es¬ 
pecially, the logarithmic negativity A f per mode is an 
increasing function of the nonclassical depth r. Also, 
the entropy S a per mode is an increasing function of the 
logarithmic participation ratio R a per mode, a = s, i. 


VII. EXPERIMENTAL MULTIMODE TWIN 
BEAMS 

Real experimental multimode twin beams have a more 
complex structure than that discussed in Sec. VI [391 
[521 S3 • The reason is that the spatiotemporal modes of 
twin beams are shared by the signal and idler fields and 
so they can be broken before or during the detection ow¬ 
ing to spectral and/or spatial filtering. As a consequence, 
real multimode twin beams are composed of three compo¬ 
nents [SUM]. A paired component describes photons em¬ 
bedded in spatiospectral modes detected by both signal- 
and idler-field detectors. A noise signal (idler) compo¬ 
nent then describes photons occurring in signal (idler) 
spatiotemporal modes that originate in filtering of the 
idler (signal) field. If we assume for simplicity that the 
paired component is ideal, i.e., without noise, we need 
six parameters to describe a real twin beam. Each com¬ 
ponent is characterized by the number M of modes and 
mean photon-pair (or noise photon) number B. The sta¬ 
tistical operator pe of the experimental twin beam can 
be expressed as 

PE = (^) Pp (^) Pn,s (^) Pn, i (37) 

M p M s Mi 

using single-mode statistical operators p p , p n<a , and p n , t 
of the photon-pair, noise signal, and noise idler compo¬ 
nents. In Eq. ( [37] ), M p , M s , and Mi give the numbers 
of equally populated modes with the mean numbers B p , 
B s , and B\ of photon pairs per mode, respectively. 

Entanglement in the experimental twin beam is cre¬ 
ated only by its paired component and as such it can be 
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quantified by the logarithmic negativity introduced 
in Sec. VI. The noise components do not contribute to 
entanglement on one side, and they do not degrade en¬ 
tanglement on the other side. This is qualitatively dif¬ 
ferent from the case of multimode twin beams discussed 
in Sec. VI and containing noise in paired spatiotemporal 
modes. 

Nonclassicality can be quantified by a multimode 
generalization of nonclassical depth te introduced in 
Ref. [IS] for a single-mode field. In a multimode twin 
beam, we may first determine the standard nonclassical 
depths T n for each single-mode field, included either in 
the paired part of the twin beam or in the noisy signal 
and idler parts of the twin beam. Then we can take either 
max„(r n ) or ^2 n r n to quantify the multimode nonclassi¬ 
cal depth te- In the first case, the nonclassical depth te 
of the experimental multimode twin beam is just given 
by the nonclassical depth r of a paired mode. The sec¬ 
ond case is physically more interesting, as the value of 
te is linearly proportional to the minimum amount of 
additional noise needed to conceal nonclassicality of the 
multimode state. In this case, we have, for the experi¬ 
mental multimode twin beams, 

t e = M p t. (38) 

Using the logarithmic negativity iV]0® defined in Sec. VI 
and the nonclassical depth te, one-to-one correspondence 
between the entanglement and the nonclassicality is ob¬ 
tained also for M-rnode twin beams. 

On the other hand, the concept of weak nonclassical¬ 
ity [ 2 U ESI m\ i s also useful for the experimental multi- 
mode twin beams considered to be composed of one effec¬ 
tive paired (macro)mode. The joint quasidistribution Pw 
of the integrated intensities W s and It 7 ,; of the signal and 
idler fields, respectively, describes the properties of this 
effective paired mode m- As no information about the 
phase is encoded in this simplified effective description, 
we may only determine the nonclassical intensity depth 
T\y quantifying nonclassicality, which demonstrates itself 
by negative values of the marginal quasidistribution of 
integrated intensities. We have to emphasize that the 
nonclassical intensity depth tw is only a nonclassicality 
witness or parameter, which reveals nonclassicality solely 
in photon-number statistics. Contrary to this, the non¬ 
classical depth r is a genuine and commonly used non¬ 
classicality measure. We note that te standard nonclas¬ 
sicality quantified by r reveals both strongly and weakly 
nonclassical states (68l [69]. From this point of view r is 
a strong tool or criterion. On the other hand, rye detects 
only strongly nonclassical states; i.e., it is a weak tool. 

The nonclassical intensity depth tw has been deter¬ 
mined for the experimental multimode twin beams in 
Ref. [M] . 

tw = Vfi 2 ~ 7 - & (39) 

where 

_ M S B S + M\B\ + 2M p Bp 
P ~ M s + Mi + 2M P ’ 


FIG. 11: (Color online) Nonclassical intensity depth tw as a 
function of the mean noise photon numbers B s and B-, for the 
mean photon-pair number B p equal to 2 [bottom, light-gray 
(yellow) area], 4 [gray (green) area], and 8 [top, dark-gray 
(blue) area], assuming M p = M s = Mi = 1. The greater the 
value of B p , the greater the value of Tw- 



7 = 


M s Bl + M ^ 2 - 2M P B P 
M s + Mi + 2M P 


(40) 


The analysis of Eq. (39) shows that the experimental 


multimode twin beam is strongly nonclassical {jw > 0 ) 
provided that 


M S R S 2 + M,Bf < 2M p B p . 


(41) 


Inequality (41) means that the multimode strong non¬ 


classicality of the twin beam is lost if the noise is suf¬ 
ficiently strong. For example, if M p = M s = Mi, 
strongly nonclassical multimode twin beams are observed 


for R 2 +B? < 2R p (see Fig. 11). This behavior is similar 


to that discussed in Sec. IV [Though the boundary given 
by Tw = 0 is quantitatively different (compare Figs. [ 6 ] 
and 11). We also have here that the greater the value of 


the mean photon-pair number B p , the greater the value 
of the nonclassical intensity depth tw- Also, the greater 
the values of mean noise photon numbers B s and R;, the 
smaller the value of the nonclassical intensity depth tw- 


Similarly as in Sec. VI, the logarithmic participation 
ratio R log can be defined for each component of the twin 
beam to quantify its dimensionality. The logarithmic 
participation ratio R log of the whole twin beam is then 
naturally given as the sum of the logarithmic participa¬ 
tion ratio R X 'm v> p of the paired component and the log¬ 


ilog 


+ 7? log 

T- -O-Mi.i 


of the noise 


arithmic participation ratio R M E 
signal and idler components. We note that Eq. (25) 
is appropriate for determining the participation ratio of 
both the single-mode noise signal (or idler) field and the 
single-mode paired field. Alternatively we may consider 
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entropies of the components instead of participation ra¬ 
tios. Entropies of the single-mode noise fields are given 
by Eq. (33). Equation (33) is applicable also for determi¬ 
nation of the entropy of entanglement of a single-mode 
paired field in a pure state for which p s jj <— c?. As a 
consequence, the entropies Sm oi a for a = p, s, i, of each 
component are increasing functions of the corresponding 
participation ratios R.M a ,a■ In single-mode cases, these 


functions are determined by Eq. (35), plotted in Fig. 10 


Similarly to the overall logarithmic participation ratio 
R log , the overall entropy S can be naturally split into 
its entangled part Sm p , p and noisy part Sm s ,s + SM it i, 
originating in the noise signal and idler components. 

Finally, we briefly address the issue of the experimental 
determination of the quantities discussed above. As these 
quantities characterize the “internal” structure of a twin 
beam, only their indirect determination is possible. It is 
based upon the measurement of the joint signal-idler pho¬ 
tocount histogram using photon-number-resolving detec¬ 
tors. Knowing these detector parameters [33], recon¬ 
struction of the joint signal-idler photon-number distri¬ 
bution El IM] provides the applied mean photon(-pair) 
numbers B and numbers M of modes. The above-derived 
formulas then give the discussed quantities. 


VIII. CONCLUSIONS 

The entanglement and nonclassicality of a single-mode 
noisy twin beam have been quantified using the negativ¬ 
ity and the nonclassical depth, respectively. Universal 
mapping between the nonclassical depth and the nega¬ 
tivity has been revealed for noisy twin beams. The map¬ 


ping reflects the fact that nonclassicality of a twin beam 
is caused by the entanglement of its two parts originat¬ 
ing in pairing of photons. Limitations to the amount 
of noise have been found to preserve entanglement to¬ 
gether with nonclassicality. the degrees of freedom of a 
twin beam quantified by the signal- and idler-field partic¬ 
ipation numbers have been divided into those needed to 
describe entanglement and the remaining ones forming 
the noisy signal and idler parts of the twin beam. The 
entanglement dimensionality derived from the negativity 
has been applied here. Entropy as an increasing function 
of the participation number has been discussed. Proper¬ 
ties of multimode twin beams have been analyzed using 
appropriate quantities related per one mode. Also, ex¬ 
perimental multimode twin beams containing additional 
noise in independent spatiotemporal modes have been in¬ 
vestigated from the point of view of their entanglement 
and multimode nonclassicality including weak nonclassi¬ 
cality and dimensionality. 


Acknowledgments 


I.A. and J.P. Jr. thank M. Bondani for discus¬ 
sions. Support by project LO1305 of the Ministry of 
Education, Youth and Sports of the Czech Republic is 
acknowledged. I.A. thanks project IGA_PrF_2015004 
of IGA UP Olomouc. J.P. Jr. acknowledges 
project P205/12/0382 of GA CR. A.M. was supported 
by the Polish National Science Centre under Grants 
No. DEC-2011/03/B/ST2/01903 and No. DEC- 
2011/02/A/ST2/00305. 


[1] A. Einstein, Ann. Phys. 17, 132 (1905). 

[2] A. Einstein, B. Podolsky, and N. Rosen, Phys. Rev. 47, 
777 (1935). 

[3] E. Schrodinger, Naturwissenschaften 23, 807 (1935). 

[4] J. Perina, Z. Hradil, and B. Jurco, Quantum Optics and 
Fundamentals of Physics (Kluwer, Dordrecht, 1994). 

[5] M. A. Nielsen and I. L. Chuang, Quantum Computation 
and Quantum Information (Cambridge University Press, 
Cambridge, 2000). 

[6] A. Migdall, Physics Today 52, 41 (1999). 

[7] V. Giovannetti, S. Lloyd, and L. Maccone, Nat. Photon. 
5, 222 (2011). 

[8] J. Perina Jr., O. Haderka, M. Hamar, and V. Michalek, 
Opt. Lett. 37, 2475 (2012). 

[9] H. Vahlbruch, A. Khailaidovski, N. Lastzka, C. Graef, 
K. Danzmann, and R. Schnabel, Class. Quantum Grav. 
27, 084027 (2010). 

[10] L. Mandel and E. Wolf, Optical Coherence and Quantum 
Optics (Cambridge University Press, Cambridge, 1995). 

[11] J. Perina, Quantum Statistics of Linear and Nonlinear 
Optical Phenomena (Kluwer, Dordrecht, 1991). 

[12] D. F. Walls and G. J. Milburn, Quantum Optics 
(Springer, Berlin, 1994). 


[13] D. Bouwmeester, J. W. Pan, K. Mattie, M. Eibl, H. We- 
infurter, and A. Zeilinger, Nature (London) 390, 575 
(1997). 

[14] G. Weihs, T. Jennewein, C. Simon, H. Weinfurter, and 
A. Zeilinger, Phys. Rev. Lett. 81, 5039 (1998). 

[15] A. Luks, V. Perinova, and J. Perina, Opt. Commun. 67, 
149 (1988). 

[16] C. T. Lee, Phys. Rev. A 42, 1608 (1990). 

[17] A. I. Lvovsky and M. G. Raymer, Rev. Mod. Phys. 81, 
299 (2009). 

[18] O. Jedrkiewicz, Y. K. Jiang, E. Brambilla, A. Gatti, 
M. Bache, L. A. Lugiato, and P. Di Trapani, Phys. Rev. 
Lett. 93, 243601 (2004). 

[19] M. Bondani, A. Allevi, G. Zambra, M. G. A. Paris, and 
A. Andreoni, Phys. Rev. A 76, 013833 (2007). 

[20] J. Perina Jr., O. Haderka, and V. Michalek, Opt. Express 
21, 19387 (2013). 

[21] M. Lamperti, A. Allevi, M. Bondani, R. Machulka, 
V. Michalek, O. Haderka, and J. Perina Jr., J. Opt. Soc. 
Am. B 31, 20 (2014). 

[22] M. Lamperti, A. Allevi, M. Bondani, R. Machulka, 
V. Michalek, O. Haderka, and J. Perina Jr., Int. J. Quant. 
Inf. 31, 1461017 (2014). 







12 


[23] R. J. Glauber, Phys. Rev. 131, 2766 (1963). 

[24] E. C. G. Sudarshan, Phys. Rev. Lett. 10, 277 (1963). 

[25] V. Dodonov and V. Manko, Theory of Nonclassical States 
of Light (Taylor & Francis, New York, 2003). 

[26] W. Vogel and D. Welsch, Quantum Optics (Wiley-VCH, 
Weinheim, 2006). 

[27] T. Richter and W. Vogel, Phys. Rev. Lett. 89, 283601 

( 2002 ). 

[28] A. Miranowicz, M. Bartkowiak, X. Wang, Y.-X. Liu, and 
F. Nori, Phys. Rev. A 82, 013824 (2010). 

[29] M. Bartkowiak, A. Miranowicz, X. Wang, Y.X. Liu, W. 
Leohski, and F. Nori, Phys. Rev. A 83, 053814 (2011). 

[30] A. Allevi, M. Lamperti, M. Bondani, J. Perina Jr., 
V. Michalek, O. Haderka, and R. Machulka, Phys. Rev. 
A 88, 063807 (2013). 

[31] W. Vogel, Phys. Rev. Lett. 100, 013605 (2008). 

[32] O. Haderka, J. Perina Jr., M. Hamar, and J. Perina, 
Phys. Rev. A 71, 033815 (2005). 

[33] J. Perina Jr., M. Hamar, V. Michalek, and O. Haderka, 
Phys. Rev. A 85, 023816 (2012). 

[34] J. Perina Jr., O. Haderka, V. Michalek, and M. Hamar, 
Phys. Rev. A 87, 022108 (2013). 

[35] C. T. Lee, Phys. Rev. A 41, 1721 (1990). 

[36] A. Verma and A. Pathak, Phys. Lett. A 374, 1009 (2010). 

[37] R. Horodecki, P. Horodecki, M. Horodecki, and K. 
Horodecki, Rev. Mod. Phys. 81, 865 (2009). 

[38] A. Peres, Phys. Rev. Lett. 77, 1413 (1996). 

[39] P. Horodecki, Phys. Lett. A 232, 333 (1997). 

[40] K. Zyczkowski, P. Horodecki, A. Sanpera, and M. Lewen- 
stein, Phys. Rev. A 58, 883 (1998). 

[41] G. Vidal and R. F. Werner, Phys. Rev. A 65, 032314 

( 2002 ). 

[42] J. F. Clauser, M. A. Horne, A. Shimony, and K. A. Holt, 
Phys. Rev. Lett. 23, 880 (1969). 

[43] E.V. Shchukin and W. Vogel, Phys. Rev. Lett. 95, 230502 
(2005); A. Miranowicz and M. Piani, Phys. Rev. Lett. 97, 
058901 (2006). 

[44] A. Miranowicz, M. Piani, P. Horodecki, and R. 
Horodecki, Phys. Rev. A 80, 052303 (2009). 

[45] C. T. Lee, Phys. Rev. A 44, R2775 (1991). 

[46] I. P. Degiovanni, M. Bondani, E. Puddu, A. Andreoni, 
and M. G. A. Paris, Phys. Rev. A 76, 062309 (2007). 

[47] M. Bondani, E. Puddu, I. P. Degiovanni, and A. An¬ 
dreoni, J. Opt. Soc. Am. B 25, 1203 (2008). 

[48] I. P. Degiovanni, M. Genovese, V. Schettini, M. Bon¬ 
dani, A. Andreoni, and M. G. A. Paris, Phys. Rev. A 79, 


063836 (2009). 

[49] C. Eltschka and J. Siewert, Phys. Rev. Lett. Ill, 100503 
(2013). 

[50] K. Zyczkowski, Phys. Rev. A 60, 3496 (1999). 

[51] M. V. Fedorov, M. A. Efremov, A. E. Kazakov, K. W. 
Chan, C. K. Law, and J. H. Eberly, Phys. Rev. A 69, 
052117 (2004). 

[52] M. V. Fedorov, M. A. Efremov, P. A. Volkov, E. V. 
Moreva, S. S. Straupe, and S. P. Kulik, Phys. Rev. Lett. 
99, 063901 (2007). 

[53] K. W. Chan, J. P. Torres, and J. H. Eberly, Phys. Rev. 
A 75, 050101 (2007). 

[54] H. Di Lorenzo Pires, C. H. Monken, and M. P. van Exter, 
Phys. Rev. A 80, 022307 (2009). 

[55] K. Bartkiewicz, K. Lemr, and A. Miranowicz, Phys. Rev. 
A 88, 052104 (2013). 

[56] F. Just, A. Cavanna, M. V. Chekhova, and G. Leuchs, 
New J. Phys. 15, 083015 (2013). 

[57] K. Bartkiewicz, J. Beran, K. Lemr, M. Norek, and A. 
Miranowicz, Phys. Rev. A 91, 022323 (2015). 

[58] K. Bartkiewicz, P. Horodecki, K. Lemr, A. Miranowicz, 
and K. Zyczkowski, Phys. Rev. A 91, 032315 (2015). 

[59] J. Perina and J. Krepelka, J. Opt. B: Quantum Semiclass. 
Opt. 7, 246 (2005). 

[60] A. Luks, J. Perina Jr., W. Leohski, and V. Perinova, 
Phys. Rev. A 85, 012321 (2012). 

[61] V. Perinova, A. Luks, J. Krepelka, and J. Perina Jr., 
Phys. Rev. A 90, 033428 (2014). 

[62] J. Pef'ina and J. Krepelka, Opt. Commun. 284, 4941 

( 2011 ). 

[63] C. K. Law and J. H. Eberly, Phys. Rev. Lett. 92, 127903 
(2004). 

[64] A. Gatti, T. Corti, E Brambilla, and D. B. Horoshko, 
Phys. Rev. A 86, 053803 (2012). 

[65] M. V. Chekhova, G.Leuchs and M. Zukowski, Opt. 
Comm. 337, 27 (2015). 

[66] D. B. Horoshko, G. Patera, A. Gatti, and M. I. Kolobov, 
Eur. Phys. J. D 66, 239 (2012). 

[67] A. M. Perez, T. S. Iskhakov, P. Sharapova, S. Lemieux, 
V. Tikhonova, M. V. Chekhova, and G. Leuchs, Opt. 
Lett. 39, 2403 (2014). 

[68] Arvind, N. Mukunda and R. Simon, Phys. Rev. A 56, 
5042 (1997). 

[69] Arvind, N. Mukunda and R. Simon, J. Phys. A: Math. 
Gen. 31, 565 (1998). 



